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O I Abstract 

A modification of perturbation theory, known as deha-expansion (vari- 
ationally improved perturbation), gave rigorously convergent series in 
Q^l some D = 1 models (oscillator energy levels) with factorially divergent 

(— I I ordinary perturbative expansions. In a generalization of variationally 

improved perturbation appropriate to renormalizable asymptotically 
free theories, we show that the large expansion orders of certain physi- 
H ' cal quantities are similarly improved, and prove the Borel convergence 

of the corresponding series for m„ < 0, with rriy the new (arbitrary) 
mass perturbation parameter. We argue that non-ambiguous esti- 
mates of quantities relevant to dynamical (chiral) symmetry breaking 
in QCD, are possible in this resummation framework. 



1 Introduction 

A "first principle" determination of tlie order parameters cliaracterizing dy- 
namical (e.g. cliiral) symmetry breaking (XSB ) in asymptotically free theo- 
ries (AFT) like QCD is traditionally considered inaccessible (except perhaps 
from lattice calculations), due to three main obstacles: 
(i) order parameters, like the quark condensate {qqY^^-, are expected to be of 
0{Kqcd)i so that the coupling at such scale is large: ordinary perturbative 
expansion is invalidated. 

(m) At arbitrary perturbative order, (gg) and other XSB order parameters 
vanish anyhow in the massless limit: their chiral limits are (perturbatively) 
trivial. 

{ill) A more subtle but equally important argument is that, attempts to 
extract genuine non-perturbative contributions to such quantities meet in- 
herent ambiguities, as indicated by the (infrared) renormalon singularities 
of perturbative expansions |]I|, 0]. Conventional wisdom thus treats {qq) and 
other non-perturbative condensates as parameters of a systematic operator 
product expansion (OPE)[], as best illustrated in the SVZ formalism|0]. 
Yet in many field theory models, definite non-perturbative results may be 
obtained from an appropriately resummed (but different) expansion, like the 
1/N expansion^, |^. There also exist powerful summation techniques, like 
the Borel method [P, |^ which, even for non Borel-summable expansions like 
in QCD typically, gives nevertheless precious informations on the nature of 
(power-like) non-perturbative contributions to a given physical quantity. An 
alternative summation method, known as delta-expansion (DE) or "varia- 
tionally improved perturbation" (VIP)0, §], is based on a reorganization of 
the interaction Lagrangian to depend on arbitrary adjustable parameters, 
to be fixed by some optimization prescription. In various models DE-VIP 
exhibits (though often rather empirically) an improved convergence of pertur- 
bative expansion. Moreover in some D = 1 models, the anharmonic oscillator 
typically, DE-VIP is equivalent to the "order-dependent mapping" (ODM) 
method [p!0[] , and optimization is equivalent to a rescaling of the adjustable 
oscillator frequency (mass) with perturbative order, which can essentially 



^Unlike the gluon condensate, the presence of XSB condensates hke {qq) in OPE's 
is however not directly inferred by infrared renormalons, these being screened by chiral 
symmetry [gj, cf. argument (ii) above. We will see that renormalons and argument (iii) 
are nevertheless relevant to XSB quantities in our context. 



suppress the factorial asymptotic behaviour of ordinary perturbative coef- 
ficients. Such a procedure was proven rigorously to converge]^, ll| (for an 



adequately rescaled mass) toward the exact result, e.g. for the oscillator 
energy levels |T^ and related quantities. 



Here we reconsider a variant of DE-VIP adapted to higher dimensional 
renormalizable theories, proposed some time ago||13||-||15||. The basic idea 



is to perform a modification of perturbative expansions in two stages: first 
exploiting specific renormalization group (RG) properties, which transform 
the ordinary expansion (in a coupling g) of certain physical quantities, de- 
pending only on g and on a mass m, in the alternative form of "mass power" 
expansions (MPE) in (m/A)° [m is the renormalization scale-invariant mass, 
A the basic RG scale and a is given by known RG coefficients]. This con- 
struction resums RG dependence to all orders (at least in specific schemes), 
and most interestingly exhibits a non-trivial massless (chiral) limit |]T^, |T5[ 
for DSB (XSB ) order parameters, or for analogous quantities like the "mass 
gap" in D = 2 modelsfl^. However, such a result turns out to be well- 
defined only in the approximation of neglecting all the purely perturbative 
(non RG) dependence. When arbitrary large orders of the complete (non- 
log) perturbative series are included, our naive mass gap result is plagued 
with ambiguities originating mainly from renormalon singularities, cf. point 
(iii) above, as we shall examine in more details here. 

However, in a second stage, an appropriate version of the (order depen- 
dently rescaled) DE-VIP can be performed on the complete MPE series in 
m/A, essentially replacing the true physical mass by an arbitrary adjustable 
mass parameter. In this note we mainly investigate the large order behaviour 
of the resulting "variational" expansion in this mass parameter^ We find that 
it produces a renormalization scheme (RS) dependent factorial damping of 
the original perturbative coefficients at large orders. Yet, unlike the oscilla- 
tor case, the damping is generally not sufficient to make the DE-VIP series 
readily convergent, in the standard perturbative regime, when the generically 
expected renormalon singularities are taken into account. But we show that 
the series can be Borel convergent, if approaching the chiral limit with the 
arbitrary mass parameter Re[rri\ < 0. These results apply formally a priori 
to any (asymptotically free) renormalizable models. Some concrete exam- 
ples are the D=2 0{N) Gross-Neveu (GN) model0 (where the mass gap is 

^See also ref. [n6| for a preliminary discussion. 



known exactly|T^); or a D = 4 gauged AFT with nj massless fermions like 
QCD, where the expected [T^ SU{nf)L x SU{nf)R — > SU{nf)v breaking is 



exhibited via non-perturbative order parameters. 



2 Transmuted mass expansion and mass gap 

In this and next section we essentially summarize some of the construction 
in ||l^-|l^. To illustrate simply the first stage, consider in a "generic" AFT 



the first RG order evolution for the renormalized "current" mass: 



70 



Ml = m{ii) [1 + 2h^g\^i) ln(Mi//i)] ^ , (1) 

where ho, 70 are one-loop RG coefficients, with 60 > for an AFT [our 

normalization is f3{g) = -h^g^ - hig^ , 7^(51) = ^y^g'^ _^ ^^^4 _^ j^ ^^^ 

the self- consistent condition Mi = m{Mi) defines Mi. Now, equivalently 
Eq. ([1|) reads 

Mi(m) = m [ln(Mi/A)]-^ = m F"^ (2) 

with A = /iexp[— l/26o5'^(/i)] the RG invariant scale, m = m{^)[2bog^{fi)]^^ 
the scale invariant mass {A = 70/(260)), and in Eq. (Q) 

F(m/A) = ln(m/A)-AlnF = AW[A'\m/AY/^] (3) 



where the Lambert |]19|| function, W[x] = Inx — Iniy, is plotted in Fig |I]. 
Eq. (0) has the remarkable property: F ~ (m/A)^/"^ for m — > 0, in contrast 
with the ordinary Log (see Fig. P, however asymptotic to F{m/A) for m ^ 
A. More precisely, on its principal branch (which is real-valued for real 
arguments), F has an alternative series expansion: 

n^)-U^r\^_^^ (4) 

of finite convergence radius Re = e~^A^. Mi(rh) in Eq. (^) thus exhibits 
different branches according to the values of the RG parameter A (see Fig. 0). 
Now, for most values of A, there is only one branch which for real rh values, 
is real and continuously matching the asymptotic perturbative behaviour of 
F at large rh: the one giving a non-zero "mass gap" Mi = A for m — » 



(region I in Fig. ^. Algebraically, the mass is obtained by expanding Eq. 
in®: 



Mi(m -^ 0) = m [{m/kf/^ + ■ ■ f 



A (1 + C(m/A) 



l/A^ 



(5) 



which may be viewed as a generalization (for m 7^ 0) of "dimensional trans- 
mutation". Note that Eq. (||) readily reproduces, e.g., the GN 0{N) model 




Figure 1: The Lambert function W{x) compared to \n.x. 

mass gap in the large A^, m — > limit (where A -^ 1 for A^ — > 00), tradi- 
tionally obtained in a different way |^. Eq. (|^) is, however, not a proof of 
dynamical XSB : rather, ^/XSB occurs, the fact that any mass is proportional 
to A is consistently incorporated by the properties of F{rh) for any rh, which 
provides an explicit bridge between the "non-perturbative" m < A regime, 
where F has power expansion (Q), and the short distance usual perturbative 
m 3> A (Log) regime. A crucial point is the difference between the usual ef- 
fective couphng g^ip^) = l/[&o ln(p^/A^)], having a Landau pole at p^ = Is?, 
and F^^{m) here, having its pole at m = 0, governing the massless limit 



(^) of the (pure RG) mass gap Eq. (0)0. Accordingly along the continuous 
branch I, Mi(m) has no singularity for < m < 00, as is clear from Eq. (^ 
and Fig. 0. 



^W{x) appears in various branches of physics ^S^ , in particular recently also in the QCD 
and RG context ||20|. Yet its connection with non-trivial chiral limit (|5|) was unnoticed 



before Il3]-[ll5|, to the best of our knowledge. 
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Figure 2: The different branches of Mi in Eq. @, for A = 4/9 (correspond- 
ing to first RG order QCD with three active quark fiavours). 



3 Pole mass gap and other DSB quantities 

Eq. (^ also defines a (lowest order) "pole" mass, being scale invariant to 
all orders (and gauge invariant as well, if gauge symmetry is relevant, as in 
QCD), thanks to its continued fraction form in Mi. Yet the genuine pole 
mass is not given simply by Eq. (|l|), as it includes non-log perturbative 
and RG contributions of arbitrary higher orders, though for most theories 
one only knowns at present the perturbative series up to the second or third 
order coefficients, like e.g. in the GN model [^, ^ or QCD[pT[]. 
A generalization of (^, perturbatively consistent [|15| with the usual pole 
mass, can be defined^: 



M^ 



[711] 



2-^ mF-^[C + F]-^ Y,dr,{2bQFy 



(6) 



n=0 



^Strictly, Eq. (|) applies only if C = &i/(26§) > 0. If C < (as in the 0{N) GN model, 
corresponding to an infrared fixed-point at g^ — — 6o/bi > 0), an alternative appropriate 
RG summation can be defined O, Ea. 



^=^>5 = 7^-A^=7Ti- (8) 



with 

F = \n[m/A]-A\nF-{B-C)\n[C + F], (7) 

26i' 26o ' 262 

F{'m) in (^, (|^) resums the RG dependence in ln[m] at two-loop order exactly 
(or even to all orders in the scheme 6j = 0, 7j = Vi > 2). Most interestingly, 
similarly to Eq. (^ F also has an (A, 5, C dependent) expansion in (m/A)^'"^ 
for sufficiently small m, with A now defined in Eq. (H). The coefficients dn 
implicitly include the non-log perturbative contributions from n-loop graphs 
(generically dominant, as discussed below), plus eventually (subdominant) 
contributions from higher RG orders. 

A similar construction can be performed for other physical quantities, at least 
those depending only on m and g. Examples are the perturbative expansion 
of the GN model vacuum energy [p!^, or in QCD the XSB order parameters 
F^/A (the pion decay constant) and (gg)(/i)/A'^|]l^, |15 . 



Now in (11), there are crucial differences with the "pure RG" mass gap, 

Eq. d): 

-The pole mass (or other physical quantities similarly) is infrared finite, 

gauge |]23[-, scale- and scheme-invariant, but the relation between the pole 

mass and e.g. the running mass is scheme dependent, which is manifested 

here by the RS-dependence in (|^) of: the perturbative coefficients dn] the 

RG coefficients A, B in Eq. (H); and of course A. 

-The dominant contributions dn in (^ behave rather generically asP^, |^: 



dn+i ^ (26o)"n! (9) 

n ^ oo 

so that the series Eq. (P) is badly divergent for any m > 0, and not even Borel 
summable: such a factorial growth of the perturbative coefficients, with no 
sign alternation, implies[0 ambiguities of 0(A), as we reexamine within the 
present context in section 5. The 0{N) GN model mass gap, at order 1/A^, 
also exhibits infrared renormalons similar p5[ to (|^), if considering only its 



naive perturbative expansion. In QCD, insertions of the (resummed) gluon 
propagator in the F.„ or {qq) perturbative expressions potentially give facto- 
rially growing asymptotic coefficients: while usually considered irrelevant in 
the m —>■ limit (cf. argument (ii) above), the factorial behaviour survives 
a priori in our construction due to the non-trivial chiral limitQ. 

^The form of those "XSB parameter renormalons" will be discussed elsewhere pSl Eq]. 



4 Variationally improved mass expansion 

We shall examine now how to possibly cure the latter potential ambiguities of 
such a resummation of DSB quantities, by combining the previous MPE series 
construction leading to e.g. Eq (^ with a specific form of delta-expansion. 
As mentioned in introduction, DE-VIP is essentially a reorganization of the 
interaction terms of the Lagrangian. More specifically here, we define a 
(linear) DE as the substitution 

m{fi) ^{1-6) m,; /(/i) ^ 6 /(/i) (10) 

within perturbative expressions at arbitrary order, where m{fi) is the renor- 
malized Lagrangian mass (in e.g. MS scheme), 6 the new expansion param- 
eter, and rriy an arbitrary adjustable mass. ([T0|) is equivalent to adding and 
subtracting to the massless Lagrangian a "trial" mass term rriy [6 interpolat- 
ing between the free {6 = 0) and the interacting massless Lagrangian {6 = 1)], 



and is entirely compatible with renormalization ||13| and gauge-invariance ||15 
The procedure then usually^ is to take the limit (5 — >■ 1 after performing a 
perturbative expansion of the relevant physical quantities to fixed order 6^, 
exhibiting a residual m^ dependence, so that an optimization prescription, 
typically the "principle of minimal sensitivity" (PMS)0, can be applied with 
respect to m^. However, we go here a step beyond this standard PMS usage 
by following more closely the logic that leads to rigorous convergence prop- 
erties of the DE method for the oscillator. 

In what follows we only investigate for simplicity the mass gap Eq (|^), but our 
construction can easily be generalized to similar DSB quantities. After ap- 
plying substitution (p!0|), M^{m,6) = J^k^^ki^)^'' can be most conveniently 
directly resummed, for 6 —>■ 1, hj contour integration[0 around 5 = 0, to 
arbitrary order K: an appropriate change of variable allowing to study the 
m{jj,) -^ (equivalently 5 — >• 1) limit in Eq. (p!0D is: 

5 = l-v/K; m^ = K"<rh^. (11) 

Eq. ([TID is simply a convenient way of parameterizing how rapidly the La- 
grangian mass ?7i(/i) -^ limit is reached (as controlled by 7 < 1) as function 
of the (maximal) delta-expansion order K. Similarly to refs. |^ the point 
is to adjust the rates at which m(/i) — ;> (5 — > 1) and K — >• 00 are simul- 
taneously reached, with no a priori need of invoking explicit optimization 



-e 



Figure 3: Singularities and equivalent integration contours in the v plane, 
for A = l. 



principle. 
The final contour integral summation takes a simple form, for K ^ oo: 

^ 1 r 
M^/A ~ V — idv e^"/™ ) F-^[v] d^ (26oFH)"" (12) 

where m = ihy/A, N is maximal perturbative order, and after deforma- 
tion the contour encircles the semi-axis Re[v] < (see Fig. ^ and also for 
simplicity we fix from now the scaling parameter in Eq. ( pT]) to its maxi- 
mal value (7 = 1) still compatible with massless limit (for ih^ — > 0). [The 
general 7 scaling ( pTj) can be analyzed |2^, ^ in a way more similar to the 
oscillator 0, [TT]], i.e. without the peculiar contour (5-summation Eq. (|I^), 
but largely complicates the algebraic analysis for renormalizable theories. In 
(|T^ we also omit some overall constant factors (due e.g. to A definition) ir- 
relevant for convergence properties, and temporarily made a RS choice such 
that B = in (P)-(||), rendering certain algebraic expressions below more 
tractable, without much loss of generality.] 

Eq. (|l^) can be well approximated analytically (at least for slightly restricted 
RS choices, as indicated above and further below): 



M^/A 



1 
2b, 



N 
'0 g=l 



N- 



E 

p=0 



T[p + q]{p + q + A){q + A) 
APT[l+p]T[l + q/A] 



p-i 



[m 



-g/A 



(13) 



where we assumed the leading renormalon behaviour Eq. (|^)0, and we used 

essentially Eq. (^) together with § dve^v^ = 27ri/r[— r] Vr. 

Now, some restrictions apply to (|I^): first, the sum over p is bounded iff 

1/A e N* (14) 

which we assume for simplicity from now. This is not much restrictive, ex- 
cept that for arbitrary AFT it is generally not possible both that A satisfies 
(|T^ and i? = in Eq. (|]), as assumed in (0). But the more general scheme 
B ^ simply makes Eq. ([T3|) algebraically more involved, without affecting 
the asymptotic behaviour and convergence properties discussed below. 
Second, strictly (|1^) is valid only asymptotically, for sufficiently large N: due 
to the finite convergence radius of expansion (^) , interchanging the sum in (^ 
and integration in ([l2|) is not rigorously justified. However, when ( |T^ ) holds, 
the formerly branch point f = is simply a pole, which allows to choose 
an equivalent contour of arbitrarily small radius around v = 0, thus always 
inside the convergence radius of (^ (see the dashed small circle contour in 
Fig. §). So, only the simple pole terms v~^ contribute to Eq. (|T2D , which 
finally sum up to (|T^). The extra contribution (around the cut at f = — e~^, 
e.g. for A = 1) gives the difference between the "exact" integral ([T2|) and 
expansion (|T3|), and can be evaluated numerically. These contributions are 
easily shown ioi A = 1 to contribute as 0{e~^^''^ ^ )h[N] relative to (p!^), 
where h[N] rapidly decreases ioi N -^ oo. Thus for large enough A^ (and/or 
small m ) those contributions are unessential for the convergence properties 
discussed below. 

The announced factorial damping of coefficients, as compared to the original 
perturbative expansion, is explicit in Eq. (|l5). Yet, closer examination indi- 
cates that the damping is insufficient to make this series for N —* oo readily 
convergent. Before considering the asymptotic behaviour of the full series 
(pTSD, it is instructive to examine the p = terms, behaving as: 

N 

^Y.n<l]/n^ + <l/A]{m")"^/^. (15) 

Q 



^The original n! coefficients in Eq. IM) correspond to r[p + q] in (13). Higher order 
refinements on infrared renormalon structure may easily be implemented: it essentially 
replaces (n— 1)! — s- r[n + 6i/(26Q)](l + ri?s/n+..) where trs depends on RS via 62 etc [||, 
without affecting the convergence properties discussed below. 



The denominators in ( pIS] ) overcompensate the numerator factorials iff 

0<A<1, (16) 

where for D > 1 AFT, A is RS dependent, as discussed in section 3. Thus, 
if our series would only consist of terms of the form Eq. (|1^), the solution 
would be simply to perform appropriate scheme changes A —>■ A in (P), (|13[) 
etc, so that a damping of coefficients larger than (or equal to) the factorial 
growth would make the series convergent. [For such RS changes in A one 
should consistently derive the corresponding change in e.g. the first few per- 
turbative coefficients di, etc, and in A, but this one-parameter RS change 
does not reintroduce any factorial behaviour in dn at large orders. Moreover, 
if (|T^) holds, any generic infrared (or ultraviolet) renormalon behaviour, of 
the form[0 ~ r" n! with r arbitrary, is damped similarly.] 
Unfortunately, the large N behaviour of (|TB|) differs from the simple "oscil- 
lator form" ([T5|), due to the p > 1 terms in expansion (H) reminiscent of 
RG properties. For any low p -^ N, renormalon factorials are still overcom- 
pensated HA < 1, but the r[l + q/A] damping decreases in strength as p 
increases, giving increasing contributions to the sum over p. All in all, the 
leading contributions to the coefficients of (p!3| ) happen at intermediate values 
of p. Nevertheless, the idea of damping factorials from appropriate RS choice 
does survive, when the series Eq. (^) is Borel transformed, as examined in 
next section. 



5 Borel convergence of DE-VIP 

A Borel integral slightly adapted to our case reads: 

BI{m) = M^{m) = 2-^m F^-^{C + F)'^ / dte-^' [1 + {2bo F)-^ V t" ] 

•^0 n=0 

(17) 

which would be (asymptotically) equal to (H) by formal expansion]] (upon 

assuming Eq. (|D), would the pole at to = 1 not make the integral (|l^ 

ill-defined. One should make a choice in deforming the contour e.g. above 

^We define the Borel transform (integrand of (O)) by dividing series coefRcients by 
(n — 1)! for convenience. Also, the summed RG-dependence rhF'^{C + F)^^ , having no 
factorial behaviour, is factored out of the Borel transform. 

10 



(or below) the pole, which results in an ambiguity, easily calculated to be 
0{e~^). Since F ~ ln[m/A] for m ^ A, an 0{A/m) ambiguity0 for the 
"short distance" (M, ?fi ^ A) pole mass is recovered. But in our construction 
Eq. (^ allows to trace the behaviour of F all the way down to m — > 0, where 
F — > 0: there the ambiguity becomes 0{1), and the naive RG-summed mass 
gap d^), which is C(A), gets an ambiguity of same order, as announced. 
Now for any given choice of contour avoiding the pole (or cut[0 at higher RG 
order) in the Borel plane t, let us apply the DE-VIP as defined in section 4, 
introducing the (5-expansion and contour resummation as in (^), this time 
on the Borel integral Eq (0). Interchanging the contour and Borel integrals, 
one can find after some algebra the asymptotic behaviour for N -^ oo: 

Jo 26o q ^[^ + q/A\ 

where we neglected for simplicity here the two-loop RG dependent {C + F)~^ 
term in ([T7|), as it does not affect asymptotic behaviour. It thus appears that 
the asymptotic behaviour of the Borel integrand in ([18|) is that of an entire 
series (at least for A > 0), i.e. with no poles for < t < oo. More pre- 
cisely, the pole at to = 1 in the original (standard) Borel integrand has been 
pushed to to — >■ +00 due to the factorial damping, so that the Borel integral 
is no longer ambiguous. However, integral (0) is not convergent, at least for 
Re[m"] > 0, so that the series is not Borel summable for standard (pertur- 
bative) m" values. 



But conversely, the integral in Eq. (18) can converge, for Re[ni'^] < 0. This 



is the case at least for A = 1, which can always be chosen by an appropriate 
and simple RS change, as previously explained. Now, since m" = rriy/A 
is an arbitrary parameter (and physical quantities anyway only depend on 
m'^ in relativistic field theories), it should be legitimate to reach the chiral 
limit m" — ;> 0, of main interest here, within the Borel-convergent half-plane 
Re[m"] < 0. For A ^ 1, one may also choose the arbitrary parameter m" 
with i?e[(m")^/'^] < such that (|18D converges, though this appears not 
possible for any arbitrary A values. This is however only an artifact of our 
simplest choice of the (5-expansion summation defining the DE-VIP series and 
leading to (|TBp: for instance, an appropriate (A-dependent) generalization of 
Eqs. (|lO])-(0) defining the DE-VIP expansion, directly leads to a Borel con- 
vergent series independently of A values |^5|. Moreover, as already mentioned 



11 



the function F{rh) in (|^) is well-defined (analytic) for any A values in a cir- 
cle of radius e~^A^ around zero (and for A = 1 the only singularity is at 
F = —1 i.e. m/A = — e^^, cf. Fig. [I|). The higher RG order F in Eq. (|^ 
has similar properties, with finite (but RS dependent |15|) convergence ra- 
dius around zero. Thus, one can choose -Re[m"] < and/or equivalently 
i?e[-F(m")] < 0, while the mass gap i?e[M(F)] = Re[M{fri)\ always remains 
positive, see Fig. ||. 

Actually, one can see directly the Borel summability of Eq. (|1^, indepen- 
dently of A, without need of considering the DE-VIP expansion: if F < 0, 
F = —\F\ simply produces the adequate sign alternation in the factorially 
growing coefficients. More precisely, a straightforward calculation of Eq. (|l^ 
for Re[F] < (again neglecting the two- loop RG dependence C, irrelevant 
to asymptotic properties), gives: 

M^/A~e-l^l + ^Ei(-|F|) (19) 

where the exponential integral function Ei{—x) has well-defined (sign alter- 
nated) asymptotic expansion for a: > 0. 

Note however that the DE-VIP expansion, leading to Eqs. (|I8[) , appears to 
improve further the asymptotic behaviour, at least for A = 1, as compared 
to (|19D, due to the extra factorial damping. For A = 1 and Re[m"] < 
Eq. (18) becomes after integration 



M^arim) ~ const. A (1 + /(|m" I) ) (20) 

where /(|m"|) — >■ exponentially fast. The (here unspecified) overall con- 
stant in ( ^0|) originates essentially from RG dependence, involving non-trivial 
factors such as 2~'^ etc at second RG order, cf. Eq.(|^), (p!7|). 
We have thus obtained Borel convergence for a certain range of the arbitrary 
mass, strictly only for Re[m"] < 0, but in which in addition the purely per- 
turbative contributions are small and even vanishing in the chiral \m"\ — > 
limit. This does not mean, though, that our final DE-VIP result is com- 
pletelly independent of the perturbative information, since the above series 
are only asymptotic to the exact series. Rather, it suggests that the "non- 
perturbative" result in the chiral limit may be essentially determined by 
pure RG properties, plus eventually the very first few perturbative terms, 
but not infiuenced by details of the large perturbative orders. Note indeed 

12 



that, only from the properties of F around -F < (thus independently of the 
DE-VIP construction) the Borel sum in (|l^) reproduces at least qualitatively 
the asymptotic behaviour of the exact 1/N result in the 0{N) GN modelQ: 
in this model the 26o in (|19D is more precisely replaced by A^ — 2, and the 
exact result [jl^ has an asymptotic expansion ||27|| similar to ([19|) , except for a 



finite term 7^;, which not suprisingly cannot be guessed by our simple Borel 
summation of the (leading) renormalon behaviour in Eq. (IT^). 

6 Discussion 

Though renormalon ambiguities are perturbative artifacts expected to dis- 
appear (or more precisely to cancel out with OPE contributions) in truly 
non-perturbative calculations [Q, ^ ^, |2^, such explicit cancellations are 
generally inaccessible for theories like QCD. Rather, the peculiar damping 
mechanism of factorial divergences exhibited here is intuitively due to the fact 
that our reorganization of perturbative expansions makes those much more 
similar to the oscillator energy levels expansion, exhibiting a dependence on 
m^j/A, Eq. (p, which is power-like (rather than log-like) for sufficiently small 
rh^. Moreover, the adjustable parameter m^/ K may be order-dependently 
rescaled, or can take arbitrary values, in particular _Re[m„/A] < producing 
sign alternation of factorial coefficients. The DE-VIP expansion appears in 
that way to "bypass" the need for explicit (and generally complicated) can- 
cellation between perturbative and non-perturbative contributions, at least 
for certain physical quantities like the mass gap. Note also that the lin- 
ear DE-VIP taking the form (p!2D, and ([I8D when combined with the Borel 



method, is only one among other possible similar resummations. In partic- 
ular, we emphasize that the obtained convergence properties for Re[F] < 
(_Re[m^,/A] < 0) do not depend on the detailed properties of the "delta- 
expansion" contour integrals here considered, e.g. Eq. ([T^p [which lead how- 
ever to a rather simple and tractable expressions in the massless limit and for 
Borel transforms Eqs. (|TBp-(^)]: more generally applying the 5-expansion 
idea in slightly different forms may replace (|13[) and subsequent results with 



eventually different series [^, but with similar asymptotic and (Borel) con- 
vergence properties. 



*The Borel summability of the exact 1/A^ 0{N) GN model mass gap, independently 
of the present construction, is analysed in details in E^ . 
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In summary, our construction exhibits an explicit counter-example to con- 
ventional wisdom arguments (i)-(ii), and to some extent (iii), mentioned in 
introduction. In the present paper we have only analyzed the formal Borel 
convergence properties, a priori applicable to any AFT, relying essentially 
on the properties of F{rh) in Eqs. (@),(^, interpolating smoothly from the 
ordinary perturbative coupling to the infrared mass power expansion. These 
convergence properties can be viewed as the generalization to D > 1 renor- 
malizable theories of the ordinary convergence properties of the DE-VIP for 
the oscillator [^ |ll|. Next we argue that such a summation recipe can pro- 
vide a well-defined basis to estimate more precisely some of the XSB order 
parameters in QCD or other models, and a more detailed study with con- 
crete numerical applications to the GN model and QCD will be explored 



in [|25[. Though one may eventually raise that in QCD-like theories, other 
contributions to the XSB order parameters of "truly non-perturbative" origin 
(i.e. unreachable by any resummation mean, and/or related e.g. to instanton 
phenomena typically) may be expected, the resummation contributions here 
considered should be a useful piece of information. 
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